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13.  ABSTRACT 


Fully  developed  turbulent  flow  of  air  at  Re  -  1.5  x  104  to  3  x  104  entered  a  porous  lube  with  circular  cross 
section.  Helium  and  air  were  injected  uniformly  through  the  tube  wall  for  36  radii  at  various  ratios  of  mass 
velocity  through  the  tube  wall  to  the  average  mass  velocity  at  the  entrance  cross  section  of  the  tube  ranging 
from  0.0001  to  0.0004.  Measurements  of  the  velocity  field  were  made  for  air  injection,  and  measurement'*  of 
helium  concentration  were  made  for  helium  injection.  Within  the  tested  range  of  air  injection,  the  distortion 
of  the  fully  developed  velocity  profile  in-terms  ot  u/u -is  negligible. 

The  numerical  solution  obdiffusion  equation  was  obtained  based  upon  the  assumption  that  the  nondimen 
sional  velocity  distribution,  u/u  ,  was  a  function  of  the  entrance  Reynolds  mtiubct  but  independent  of  small 
wall-mass  injection  rates.  Prandtl  mixing  length  theory  incorporating  Van  Driest’s  damping  factor  generalised 
by  Kinney  and  Sparrow  for  the  case  of  wall-mass  transfer  was  used  to  evaluate  effective  diffusivity  for 
turbulent  flow.  The  accuracy  of  the  present  theory  was  compared  with  Kays’  correlation  of  turbulent  heat 
transfer  data  with  zero  mass  injection,  Raithby  ’s  exact  solution  of  laminar  heat  transfer  with  and  without 
wall-mass  transfer,  and  the  present  experiments  of  turbulent  diffusion.  Excellent  agreement  with  present 
theory  is  obtained  in  all  these  comparisons. 

For  convenience  of  users  of  the  present  theory ,  complete  computer  listing  w ith  users’  instruction  is  appended 
to  the  report.  The  results  arc  pertinent  to  reverse  osmosis  desalination  processes,  water  supply  systems, 
catalytic  converters  for  pollutant  control,  and  the  numerous  subject  areas  in  which  a  mixing  of  gases  is 
predominant. 


14 


UNCLASSIFIED 

Stcunty  Classification 


KEY  WOROl 


Diffusion 

Heat  and  Mass  Transfer 
Evaporation 
Absorption 
Simulation 

Turbulent  Boundary  Layer 
Laminar  Boundar)  Layer 
Injection  of  Cases 
Computer  Programs 


U.  S.  ARMY  MOBILITY  EQUIPMENT 
RESEARCH  AND  DEVELOPMENT  CENTER 
FORT  BELVOIR,  VIRGINIA 


Report  2053 


TURBULENT  DIFFUSION  IN  ROUND  TUBES 

Task  1T06J 101A91A00 


February  1973 


Distributed  by 
The  Commander 

U  S.  Army  Mobility  Equipment  Research  and  Development  Center 


Prepared  by 

J.  A.  Zwolinski 
S.  W.  Chi 
A.  S.  Pa  til 

Independent  Laboratory  In-House  Research 
Mechanical  Technology  Department 


Approved  for  public  release;  distribution  unlimited. 

iJ> 


SUMMARY 


Experiments  were  performed  with  fully  developed  turbulent  flows  of  air  at  Re  = 
1.5  x  104  to  3  x  104  entering  a  porous  tube  with  circular  cross  section.  Helium  and  air 
were  injected  uniformly  through  the  tube  wall  for  36  radii  at  various  ratios  of  mass  ve¬ 
locity  through  the  tube  wall  to  the  average  mass  velocity  at  the  entrance  cross  section 
of  the  tube  ranging  from  0.0001  to  0.0004.  Measurements  of  velocity  field  were  made 
for  the  case  of  air  injection,  and  measurements  of  helium  concentration  were  made  for 
the  case  of  helium  injection.  Within  the  tested  range  of  air  injection,  the  distortion  of 
the  fully  developed  velocity  profile  in  terms  of  u/u  is  negligible. 

Numerical  solution  of  the  diffusion  equation  was  obtained  based  upon  the  assump¬ 
tion  that  the  nondimensional  velocity  distribution,  u/u,  was  a  function  of  the  entrance 
Reynolds  number  but  independent  of  small  wall-mass  injection  rates.  Prandtl  mixing 
length  theory  incorporating  Van  Driest’s  damping  factor  generalized  by  Kinney  and 
Sparrow  for  wall-mass  transfer  was  used  to  evaluate  effective  diffusivity  for  turbulent 
flow.  Accuracy  of  the  present  theory  was  compared  with  Kays’  correlation  of  turbu¬ 
lent  heat  transfer  data  with  zero  mass  injection,  Raithby’s  exact  solution  of  laminar 
heat  transfer  with  and  without  wall-mass  transfer,  and  the  present  experiments  of  tur¬ 
bulent  diffusion.  Excellent  agreement  with  present  theory  was  obtained  in  all  these 
comparisons. 

For  the  convenience  of  users  of  the  present  theory,  a  complete  computer  listing 
with  users’  instruction  is  appended  to  the  report.  The  results  are  pertinent  to  reverse 
osmosis  desalination  processes,  water  supply  systems,  catalytic  converters  for  pollutant 
control,  and  the  numerous  subject  areas  in  which  a  mixing  of  gases  is  predominant. 
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Equation  (21) 
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TURBULENT  DIFFUSION  IN  ROUND  TUBES 


I.  INTRODUCTION 


1.  The  Problem  Investigated.  Figure  1  illustrates  a  porous  tube  with  a  main¬ 
stream  of  gas  flowing  in  the  axial  direction.  Another  gas  may  flow  through  the  pores 
in  the  tube  surface  to  join  the  mainstream.  In  addition,  liquid  droplets  may  be  injected 
into  the  gas  stream.  The  flow  in  the  tube  is  turbulent.  Large  property  variations  may 
exist  in  the  fluid  as  a  result  of  one  or  more  of  the  following  influences:  (1)  temperature 
differences  existing  between  the  surface,  the  mainstream,  the  secondary  stream,  and  the 
liquid  droplets;  (2)  composition  differences  between  the  different  fluids;  and  (3)  chem¬ 
ical  reactions  between  one  or  more  of  the  fluids. 


liquid 
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Fig.  1.  Diagrammatic  representation  of  the  general  problem  of  heat  and  mass 
transfer  in  round  tubes. 


The  flow  represented  in  Fig.  1.  contains  features  which  are  present,  though 
not  usually  simultaneously,  in  a  great  many  circumstances  of  practical  importance.  If 
only  one  fluid  is  present  (no  injection  of  liquid  droplets  or  secondary  gas),  the  situation 
is  the  familiar  one  of  turbulent  tube  flow;  it  arises  in  heat  exchangers  and  numerous 
other  engineering  devices.  Injection  or  suction  of  fluid  through  the  porous  tube  wall 
occurs  in  transpiration  cooling  of  the  wall  of  gas-turbine  combustors  and  reverse-osmosis 
desalination  of  saline  water;  injection  may  also  be  resorted  to  as  a  means  to  prevent 
boundary-layer  separation  (e.g.,  to  improve  efficiency  of  compressor  diffusers).  A  chem 
ical  reaction  exists  between  the  secondary  and  mainstreams  in  automobile  afterburners 
and  catalytic  converters  for  pollutant  control.  Interaction  between  liquid  droplets  and 
the  gas  stream  occurs  in  such  practical  processes  as  evaporation  of  water  droplets  in 
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humidifiers,  absorption  of  toxic  gas  by  liquid  droplets  in  environmental  control  equip¬ 
ments,  and  combustion  of  fuel  droplets  in  various  thermal  power  units. 

The  ultimate  aims  of  the  present  program  at  USAMERDC  are  the  provision 
of  a  theory  for  the  prediction  of  heat,  mass,  and  momentum  transler  in  situations  such 
as  those  of  Fig.  1  and  applications  of  the  theory  to  optimize  performance  of  a  variety 
of  mobility  equipments.  However,  the  first-year  investigations  reported  herein  have 
been  concentrated  mainly  on  turbulent  diffusion  in  round  lubes.  These  results  are  per¬ 
tinent  to  automobile  afterburners,  catalytic  converters,  reverse  osmosis  desalination  sys¬ 
tems,  and  water-supply  systems. 

2.  Review  of  Present  State  of  Knowledge.  Whereas  numerous  analytical  solu¬ 
tions  have  been  obtained  for  Laminar  diffusion, ,'s  limited  analyses  have  been  made  for 
turbulent  diffusion  in  tubes.  These  included  an  analysis  of  salt  diffusion  in  water  by  the 
method  of  analogy  between  momentum  and  mass  transfer,6  and  an  analysis  of  the  fully 
developed  concentration  profiles.7  Analyses  for  the  development  of  concentration  pro¬ 
files  for  turbulent  tube  flow  have  not  yet  been  made. 

Tests  for  turbulent  velocity  and  temperature  distributions  have  been  made 
by  Weissberg  and  Berman8,  Yuan  and  Galowin9,  Olson  and  Eckert10,  and  many  others. 

A  good  review  is  found  in  Olson  and  Eckert’s  paper11.  No  known  experimental  results 
for  concentration  distribution  in  tubes  have  been  reported. 


*G.  Raithby,  “Laminar  Heat  Transfer  in  the  Thermal  Entrance  Region  of  Circular  Tubes  and  Two-Dimensional 
Rectangular  Ducts  with  Wall  Suction  and  injection,”  Int.  J.  Heat  Mass  Transfer,  Vol.  14,  pp.  223-243  (1971). 

2 

U.  Merten,  11.  K.  Lonsdale,  and  R.  L.  Riley,  “Boundary  Layer  Effects  in  Reverse  Osmosis,”  /nd.  Eng.  Chem. 
Fundamentals,  Vol.  3,  No.  3,  pp.  210-213  (1964). 

^T.  K.  Sherwood,  P.  L.  T.  Brian,  R.  E.  Fisher,  and  L.  Dresner,  “Salt  Concentration  at  Phase  Boundaries  in  Desali¬ 
nation  by  Reverse  Osmosis,”  Ind.  Eng.  Chem.  Fundamentals,  Vol.  4,  No.  2,  pp.  113-118  (1965). 

%.  N.  Gill,  C.  Ticu,  and  D.  W.  Zch,  “Concentration  Polarization  Effects  in  a  Reverse  Osmosis  System,”  Ind.  Eng. 
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9T.  K.  Sherwood,  P.  L.  T.  Brian,  R.  E.  Fisher,  and  L.  Dresner,  “Salt  Concentration  at  Phase  Boundaries  in  Desali¬ 
nation  by  Reverse  Osmosis,  Ind.  Eng.  Chem.  Fundamentals,  Vol.  4,  No.  2,  pp.  113-118  (1965). 
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R.  B.  Kinney  and  E.  M.  Sparrow,  “Turbulent  Flow,  Heat  Transfer,  and  Mass  Transfer  in  a  Tube  with  Surface 
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It  may  be  concluded  from  the  above  review  that  knowledge  of  turbulent  dif¬ 
fusion  in  tubes  had  been  defective  because  drastic  assumptions  are  made  in  theoretical 
analyses  and  no  experiments  are  available  to  verify  the  theory. 

3.  Purpose  and  Scope  of  Present  Work.  This  work  is  intended  to  remedy  the 
above  defects.  For  this  purpose,  experimental  data  have  been  obtained  for  turbulent 
diffusion  of  helium  in  a  round  tube,  a  theoretical  model  has  been  developed  for  calcu¬ 
lation  of  turbulent  diffusions  in  round  tubes,  and  confidence  in  the  present  theoretical 
model  has  been  established  by  comparison  with  the  experiments. 

Descriptions  of  the  experiments  and  theoretical  development  are  given  in 
Sections  II  and  III.  Section  IV  is  devoted  to  a  discussion  of  the  results.  Conclusions 
for  applications  of  the  present  theory  aie  described  in  Section  V.  For  the  convenience 
of  the  users  of  the  present  theory,  a  complete  listing  of  the  computer  program  for  the 
present  theory  is  appended  to  this  report. 


II.  EXPERIMENTS 


4.  Apparatus  and  Instrumentation.  The  apparatus  used  for  the  xperiments 
consisted  of  a  flow  circuit,  an  air  blower,  a  hot  film  anemometer  system,  a  gas  analyzer 
system,  a  traversing  mechanism,  and  a  metered  helium  supply  (Fig.  2). 


POTENTIOMETER  MANOMETER  PITOTSTATIC  FLOWMETER  SAFETY  LIGHT  ROTOMETER 


HOT  FILM  100  CFM  ENTRANCE  DUCT  TEST  SECTION  WITH  EXIT  DUCT  TRAVERSING  GAS  ANALYZER 

ANEMOMETER  BLOWER  SHROUD  AND  INLET  MECHANISM 


PORTS  FOR  SECOND¬ 
ARY  GAS 
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Fig.  2.  Layout  of  the  diffusion  apparatus. 
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a.  Flow  Circuit  and  Blower.  A  schematic  diagram  of  the  flow  system  is 
shown  in  Fig.  3.  Air  was  supplied  from  a  pressure  blower  with  its  flow  rate  controlled 
by  a  gate  valve  upstream  of  the  blower.  Downstream  of  the  blower  is  an  entrance  tube 
1.875  inches  ID  and  6  feet  long  followed  by  a  porous  test  tube  1.875  inches  ID  and  3 
feet  long.  A  traversing  mechanism  was  mounted  at  the  downstream  end  of  the  test 
tube;  a  description  of  the  traversing  mechanism  is  given  in  paragraph  4d. 

Secondary  gas  helium,  metered  by  an  SK  Rotameter,  was  first  intro¬ 
duced  into  a  shroud  which  surrounds  the  test  tube  and  then  injected  into  the  main¬ 
stream  air  through  the  porous  wall  of  the  test  tube.  The  shroud  was  fi.'ed  loosely  with 
fiberglass  to  insure  pressure  uniformity.  Fine  grade  (5-ja  pore  size)  sintered  stainless 
steel  test  tube  supplied  by  Pall  Trincor  Corporation  was  used  for  wall  smoothness.  Be¬ 
cause  of  the  longitudinal  pressure  drop  inside  the  porous  tube  with  axial  flow  through 
it,  the  radial  pressure  drop  through  the  tube  wall  itself  should  be  large  so  that  the  radial 
injection  will  be  uniform  along  its  length.  Sintered  stainless  steel  tube  is  too  permeable 
to  insure  this;  20  layers  of  rayon  cloth  were  wrapped  around  the  porous  tube.  With 
this  combination,  a  pressure  differential  of  about  2  psi  is  produced  at  a  radial  helium 
flow  rate  of  1  SCFM/sq  ft. 

The  above  described  flow  system  enaLled  tests  to  be  conducted  with  a 
fully  developed  turbulent  flow  of  air  in  a  smooth  tube  entering  into  a  porous  tube  (36 
diameter  long)  at  Re  =  1.5  x  104  to  3  x  104  with  uniform  helium  injection  along  the 
porous  tube. 

Vertical  traversing 
mechanism 


Fig.  3.  Schematic  representation  of  the  flow  system. 
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b.  Hot-Film  Anemometer  System.  A  hot-film  anemometer  system  pur¬ 
chased  from  Lintronic  Laboratories  was  used  to  measure  the  velocity  distributions  for 
turbulent  flow.  The  system  consisted  of  a  standard  wedge  probe  and  a  linear  hot-film 
anemometer  with  its  output  monitored  by  a  10-v  d.c.  digital  voltmeter. 

The  probe  stem  was  2  mm  in  diameter.  The  sensing  element  was  1  mm 
long  and  0.2  mm  across,  and  the  wedge  angle  was  8  degrees.  The  length  of  the  sensing 
element  was  set  normal  to  the  tube  axis.  The  probe  was  connected  to  the  Lintronic 
Constant  Temperature  Anemometer  (Model  40B). 

The  air  velocity  and  system  output  voltage  relation  was  linear.  A  typi¬ 
cal  linear  calibration  curve  is  shown  in  Fig.  4.  The  calibration  was  performed  by  a 
standard  Aero-Lab  pitot-static  tube  and  an  NIL  precision  manometer  reading  to  the 
nearest  0.001  inch. 


Fig.  4.  Calibration  of  the  Lintronic  anemometer. 
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c.  Gas  Analyzer  System.  A  GOW-M AC  Model  20-150  thermal  conductiv¬ 
ity  Analyzer,  with  a  diffusion-type  thermal  conductivity  cell,  was  used  to  measure  he¬ 
lium  concentration.  The  cell  contained  a  sample  and  reference  gas  flow  geometry.  Two 
detector  elements,  which  formed  an  electrical  Wheatstone  bridge  circuit  and  were  heat¬ 
ed  by  an  electrical  current  from  a  constant  voltage  d.c.  power  source,  were  installed  in 
each  flow  system. 

The  instrument  was  balanced  by  passing  reference  gas  (air)  through  both 
the  sample  and  reference  flow  systems.  When  a  sample  gas  of  different  composition 
was  introduced  into  the  sample  system,  a  change  in  the  rate  of  heat  loss  occurred  caus 
ing  the  Wheatstone  bridge  to  unbalance.  The  output  voltage  from  the  detector  cell  was 
fed  into  an  analog  readout.  The  instrument  was  calibrated  by  a  standard  helium/air 
mixture.  Figure  5  shows  a  typical  calibration  of  helium  concentration  versus  analog 
readout. 


0  20  40  60  80  100 

Analog  readout 


Fig.  5.  Calibration  of  the  GOW-MAC  gas  analyzer  for  helium  in  air. 
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d.  Traversing  Mechanism.  Figure  6  shows  the  traversing  mechanism  on 
which  a  gas  sampling  probe  was  mounted.  The  gas  sampling  probe  was  ieplaccd  by  a 
hot-film  probe  when  velocit)  measurements  were  made.  Figure  7  shows  the  hot-film 
probe  laid  between  the  traversing  mechanism  and  the  Lintronic  anemometer. 


V5219 

Fig.  6.  Traversing  mechanism  with  gas  analyzer. 
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Fig.  7.  Traversing  mechanism  with  anemometer. 
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The  traversing  mechanism  allowed  probes  to  travel  in  both  the  axial 
and  the  radial  directions.  The  axial  distance  of  the  probe  tip  from  the  upstream  end  of 
the  porous  tube  was  indicated  by  a  Cenco  optical  bench.  The  radial  position  of  the 
probe  relative  to  the  tube  wall  was  indicated  by  a  vernier  slider  reading  to  the  nearest 
0.005  cm.  Closing  of  an  electric  circuit,  indicated  by  lighting-up  of  a  torch  bulb,  was 
used  to  signal  the  contact  of  the  probe  with  the  tube  wall. 

5.  Test  Procedure  and  Measurements. 

a.  Velocity  Measurements  With  and  Without  Wall-Mass  Injection.  The 
main  air  flow  was  measured  by  a  pitot-static  tube  placed  centrally  at  x/D  =  25  from 
the  leading  edge  of  the  entrance  tube.  The  pressure  drop  and  the  static  pressure  were 
measured  with  an  NIL  precision  manometer  reading  to  0.001  in.  of  water.  The  temper¬ 
ature  at  the  exit  was  measured  by  a  thermometer  reading  to  0.2°  F.  From  tlwsc  data, 
centerline  velocity  of  the  main  air  flow  was  obtained.  The  axial  velocity  u,  was  mea¬ 
sured  by  a  hot  film  anemometer  at  10  essentially  equal  intervals  of  (r/R)1  at  (x/D)  =  0, 

5,  10,  and  20  from  the  leading  edge  of  the  porous  tube. 

The  above  test  runs  were  made  at  main  flow  Re  =  1.5  x  104.  Runs 
were  made  first  without  mass  injec  tion  at  the  porous  tube  and  then,  with  uniform  air 
injection  along  the  porous  tube  at  m  (=  pvw/pu0)  =  -0.0004. 

b.  Concentration  Measurements  with  Helium  Injection.  The  concentration 
of  helium  distribution  was  measured  with  a  gas  analyzer.  Tests  were  conducted  with  a 
fully  developed  turbulent  flow  of  air  in  the  entrance  tube  entering  the  porous  tube  at 
Re  =  1.5  x  104  and  3  x  104  with  uniform  helium  injection  along  the  porous  tube  at 

m  (=m"nc/pu0)  =  -0.0001,  -0.0002,  and  -0.0004. 

For  each  run,  helium  concentration,  p{Ic/p,  was  measured  at  x/D  (from 
leading  edge  of  the  porous  tube)  =  5, 10,  and  20  and  at  about  15  different  r/R  positions. 

6.  Test  Results. 

a.  Velocity  Profiles.  The  measured  velocity  profiles  for  the  case  of  no 
mass  injection  at  the  porous  tube  are  plotted  in  Fig.  8  in  the  form  of  u/u0  versus 
i?(=r/R).  The  measured  velocity  profile  with  air  injection  at  the  porous  tube  is  also 
plotted  in  Fig.  8  in  the  form  of  u/u  versus  17. 

b.  Helium  Concentration  Profiles.  The  measured  helium  concentrations, 

f  (=pJ(c/p),  are  plotted  versus  dimensionless  radius,  r ?,  in  Figs.  9,  10,  and  11  for  different 
values  of  Re  (from  1.5  x  104  to  3  x  I04),  x/D  (from  0  to  20),  and  m  (from  -0.0001  to 
-0.0004). 
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0  0.2  0.4  0.6  0.8  1.0 

v  (=  r/R) 

Fig.  8.  Measured  turbulent  velocity  profiles  at  Re  =  1.5  x  104 :  X,  m  =  0;  0,  m  =  -0.0004. 
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lelium  concentration  at  He  =  1.5  x  104  and  m  =  -0.0002: 
'D  =  10;  A,  x/D  =  5. 


Helium  concentration 


III.  THEORETICAL  DEVELOPMENT 


7.  Theoretical  Model. 

a.  Diffusion  Equation.  The  problem  analyzed  is  the  diffusion  of  the  in¬ 
jected  gas,  helium  for  the  present  investigation,  in  the  fully  developed  turbulent  air 
flowing  in  a  circular  tube  as  represented  by  Fig.  12.  The  conservation  equation  for  the 
injected  helium  can  be  written: 

pu  3f/3x  +  pv  3f/3r  =  (1/r)  3(rpDc  3f/3r)  /  3r .  (1) 

This  equation  considers  convection  in  the  axial  direction  and  diffusion  and  convection 
in  the  radial  direction  but  neglects  axial  diffusion.  The  normalized  concentration,  f,  is 
the  local  helium  density  divided  by  the  mixture  density.  The  latter,  following  the  usual 
assumption  of  incompressibility,  is  considered  constant.  The  boundary  conditions  to 
be  imposed  upon  the  solution  of  Equation  (1)  are 


At  x=0,  any  r;  f=0 

(2) 

Atr=0,  anyx:  3f/3r=0 

(3) 

At  r=R,  any  x:  3f/3r=(l-f)m"w/pD 

(4) 

Fig.  12.  Nomenclature  for  flow  through  a  porous  tube  with  injection. 
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Equation  (2)  represents  the  condition  that  the  helium  concentration  is  uniformly  zero 
at  the  entrance  to  the  leading  end  of  the  porous  tube;  Equation  (3)  implies  axisymme- 
try  with  respect  to  the  tube  centerline;  and  Equation  (4)  relates  helium  diffusion  and 
convection  at  the  porous  tube  surface  to  helium  injection  rate  per  unit  area,  n\v  •  ^  *s 
also  implied  in  Equation  (4)  that  tui  julence  is  completely  damped  at  the  tube  surface; 
thus,  the  effective  diffusivity ,  De,  at  r=R  is  replaced  by  the  molecular  Uiffusivity  of  he¬ 
lium  in  air,  D. 


To  integrate  Equation  (1)  with  boundary  conditions,  Equations  (2) 
through  (4),  the  turbulent  velocity  field,  u  and  v,  and  the  effective  diffusivity,  De,  must 
first  be  specified.  Consideration  is  given  below  to  the  mathematical  representation  of 
these  quantities. 

b.  Turbulent  Velocity  Field.  A  complete  solution  for  the  velocity  field  re¬ 
quires  numerical  integration  of  continuity  and  momentum  equations.  However,  when 
the  nondimensional  mass  injection  rate,  m(=mw"/pu ),  is  small,  it  can  be  expected  that 
the  fully  developed  turbulent  velocity  profile  at  the  entrance  to  the  porous  tube  will 
not  be  distorted.  It  seems  reasonable  for  the  present  investigation  at  m  =  -0.0001  to 
-0.0004  to  assume  that  the  axial  velocity  profile  can  be.  represented  by 

u  =  u0o(r?)  (5) 

where  0O  (7?)  is  the  fully  developed  turbulent  velocity  profile,  u/u  ,  at  entrance  to  the 
porous  tube,  i.e.  at  x=0,  and  u  at  the  other  axial  position  can  be  calculated  by  the 
equation 


u=u0  (l+2m£) .  (6) 

Substitution  of  Equations  (5)  and  (6)  in  the  continuity  equation 

3ur/3x+3vr/3r  =  0  (7) 

followed  by  integration  yields  the  following  expression  for  the  radial  velocity  distribution 


v 


(-2mu/i?) 


0O  Vdv  ■ 


(8) 


Hence  u  and  v  distribution  for  specified  uQ  and  m  can  be  evaluated  if  0o(i 7)  is  determined. 

It  is  shown  in  Appendix  A  that  0o(i?)  can  be  evaluated  by  the  equation 
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0o  (tj)  =  |  J1  2R+tj  |l+  1 1+4;?  (R+)2  |  (0.14-0.08t?2-0.06t?4)  (1-exp  (-R+ 
(i-t?)/a))|2  |,/2|  2R+t?  |l+  |  l+4r?  (R+)2  |  (0. 
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-0.08t?2-0.06t74)  (l-exp(-R+(l-T?)/A))  r  n  dr?  dr? 


(9) 


where  R+  is  related  to  the  entrance  Reynolds  number.  Re.  This  relation  is  shown  graph¬ 
ically  by  Fig.  13.  Figure  14  shows  the  calculated  <j>0  (r?)  at  Re  =  2  x  103  to  l  x  106. 

c.  Effective  Diffusivity  for  Turbulent  Flow.  The  effective  diffusivity,  Dc, 
can  be  obtained  by  use  of  a  formula12  of  the  form 


Dc  =  p/pSc+pt/pSct 


(10) 


where  p  is  the  density  of  the  air,  p  is  the  laminar  viscosity  of  the  air,  Sc  is  the  laminar 
Schmidt  number  of  helium  in  air,  Sct  is  the  turbulent  Schmidt  number  which  is  a  con¬ 
stant  of  0*86,  and  pt  is  the  turbulent  viscosity.  Among  the  various  candidate  represen¬ 
tations  for  the  turbulent  viscosity,  a  model  involving  PrandtPs  mixing  length  and  Van 
Driest’s  damping  factor  has  gained  broad  acceptance  in  recent  years.13"18 

This  approach  is  employed  here  and  gives  the  following  model  for  the 
turbulent  viscosity: 


12S.  V.  Patankar  and  D.  B.  Spalding,  “Heat  and  Mass  Transfer  in  Boundary  Layers,”  Morgan  and  Grampin,  London 
(1967). 

^R.  B.  Kinney  and  E.  M.  Sparrow,  “Turbulent  Flow  Heat  Transfer,  and  Mass  Transfer  in  a  Tube  with  Surface 
Suction,”  ASME  Paper  No.  69-HT4(1969). 

*4E.  M.  Sparrow,  V.  K.  Jonsson,  G.  S.  Beavers,  and  R.  G.  Over,  “Incompressible  Turbulent  Flow  in  a  Permeable 
Wall  Duct,”  ASME  Paper  No.  71-WA/FE-l  (1971). 

^T.  Ccbcci,  A.  M.  0.  Smith,  and  G.  Mosinskis,  “Solution  of  the  Incompressible  Turbulent  Boundary-Layer 
Equations,”/.  Heat  Transfer,  Vol.  92,  pp.  133-143(1970). 

16S.  V.  Patankar  and  D.  B.  Spalding,  “Heat  Transfer  in  Boundary  Layers,”  Morgan  and  Grampin,  London  (1967). 

17S.  W.  Chi  and  C.  C.  Chang,  “Effective  Viscosity  in  a  Turbulent  Boundary  Layer,”  AlAA  Journal,  Vol.  7,  pp. 
2032-2035  (1969). 

*®W.  j.  Glowacki  and  S.  W.  Chi,  "Effect  of  Pressure  Gradient  on  Mixing  Length  for  Equilibrium  Turbulent 
Boundary  Layers,”  AIA A  Paper  No.  72-213(1972). 
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Fig.  13.  Reynolds  number,  Re,  versus  nondimensionai  tube  radius,  R+. 

=  pl2/3u/3r/  (11) 

1/R  =  DF(0.14-0.08t?2-0.06t?4)  (12) 

DF  =  1-exp  j-x-lA/2  [  (x4+4  (R+  (l-rj)/A)4]’/2  +  X2  }  **  (13) 

where  x  is  defined  as  mRe  (l-rj)/4.  Equation  (11)  expresses  Prandtl’s  mixing  length 
concept,  while  Equation  (12)  is  a  representation  for  the  mixing  length,  1,  deduced  by 
Nikuradse19  from  measured  velocity  profiles  for  the  circular  tube.  Van  Driest’s  damp¬ 
ing  factor20  generalized  by  Kinney  and  Sparrow21  is  expressed  by  Equation  (13). 


^11.  Schlichting,  “Boundary  Layer  Theory,”  McGraw-Hill,  New  York  (1960). 

2®E.  R.  Van  Driest,  "On  Turbulent  Flow  Near  a  Wall,”/.  Aero.  ScL,  Vol.  23,  pp.  1006-1011, 1036(1956). 

21r.  b.  Kinney  and  E.  M.  Sparrow,  “Turbulent  Flow,  Heat  Transfer,  and  Mass  Transfer  in  a  Tube  with  Surface 
Suction,"  ASME  Paper  No.  69-HT-4(1969.). 
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Combining  Equations  (10)  through  (13),  followed  by  rearrangement, 
yields  a  De  (=De/t>)  expression: 

De  =  l/Co+  [0.5  Re  (l+2m£)/Sct]  j  [1-exp  |-X-lA/2  [  (x4+4  (R+  (1-7?) 

m )4  ] /z+X2  (0.14-0.08t?2  -0.06r?4)  | 2  •  I  d0o/dr?l  (14) 

d.  Nondimensional  Diffusion  Equation.  The  diffusion  Equation  (1)  and 
its  boundary  conditions  represented  by  Equations  (2)  through  (4)  can  now  be  trans¬ 
formed  into  the  following  nondimensional  form: 


U(9f/9|)  +  V(9f/9i?)  =  (1/r?)  9  (7?De9f/97?)/97? 

(15) 

At  £=0,  any  7?:  f=0 

(16) 

At  77=0,  any  £:  9f/97?=0 

(17) 

At  77=1 ,  any  £:  9f/97?=0.5  Rem  Sc(l-f) 

(18) 

where  U=0.5Re  (l+2m£)  0O 

(19) 

rv 

V=  (-Re  m/7?)  I  7?0O  d7? 

J  0 

(20) 

and  Dc  and  0O  are  evaluated  by  Equations  (9)  and  (14)  respectively. 

8.  Solution  of  Diffusion  Equation  by  Finite  Difference. 

a.  Finite  Difference  Scheme.  A  two-layer  implicit  finite  difference  proce¬ 
dure  given  by  Crank  and  Nicholson22  is  developed  to  solve  the  dimensionless  diffusion 
Equation  (15)  with  boundary  conditions  of  Equations  (16)  to  (18).  Values  of  U,  V, 
and  Dc  from  known  values  of  Sc  and  m  are  first  calculated  by  Equations  (19),  (20), 
and  (14)  respectively.  Using  a  rectangular  grid  system  in  (£,  7?)  plane  with 


*i+l  =  *i  +  ^ i 

(21) 

7?j+l  =  l?j  +  At?j  , 

(22) 

00 

J.  Crank  and  P.  Nicholson,  “A  Practical  Method  for  Numerical  Evaluation  of  Solutions  of  Partial  Differential 
Equations  of  Heat  Conduction  Type,”  hoc.  Cambridge  Phil  Soc.,  Vol.  43,  p.  50  (1967). 
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partial  differential  Equation  (15)  is  approximated  at  the  midpoint  (£i+J /2.  i?j)  as  follows: 

«il'H+^I+Vj*l=si  <23> 

where 


=  De)j.jT  /  2*?j  ( At?j)2 -V-  Vi  /  4(Arjj) 

/3j  =  Uj+V2  /  +  [(7j  De):;^  +  (r;  D c)^  ]  /  2^  (A^)2 

V2  —  i+/4 

7i  =  Vj  /4(Ar?j)-(r?De)j+!/2/2I3(A7?j)2  . 

i  __  S  +  V?  1+ ^ 

=  fj_i  [  (*?  Dc)j.j,-.  /  2i?j  (Ar?j)2  +  /  4  (Arjj)  ] 


«  /  ; f >4  —  i+54  — 

+fj  |Uj  /  AS;  -  [  (77  De)j+i/2  +  (t?  Dc)j_i/2  ]  /  2r?j  (AT7p  j 

<27> 

The  grid  point  spacings  are  written  as  A£j  and  At?-  to  indicate  that  they 
can  be  varied  in  the  different  parts  of  the  grid.  Coefficients,  Oj,  0-,  y-,  and  8-  written 
above  are  valid  for  Arjj_j  =  At?j.  In  the  present  application,  however,  it  was  found  best 
to  use  a  small  interval  near  the  wall  and  to  halve  the  interval  a  number  of  times  as  77  in¬ 
creased  as  shown  in  the  Table. 

Variation  of  At? 


j 

V 

A77 

0-  39 

0  to  0.512 

0.0128 

40-  79 

0.512  to  0.768 

0.006'4 

80-119 

0.768  to  0.896 

0.0032 

120-159 

0.896  to  0.960 

0.0016 

160-199 

0.960  to  0.992 

0.0008 

200-209 

0.992  to  0.996 

0.0004 

210-219 

0.996  to  0.998 

0.0002 

220-240 

0.998  to  1.000 

0.0001 

If  the  interval  is  halved  at  successive  values  of  j,  i.e.,  Arjj_j  =  2  Arjj  = 

Arjj  +  Ar?j  +  2 ,  the  appropriate  expressions  for  the  coefficients  ctj,  ft,  7j,  and  5j  are  ob¬ 
tained  from  Equations  (24)  through  (27)  by  replacing  j+1  with  j+2. 

The  computing  procedure  involves  calculation  of  the  concentration  pro¬ 
files.  At  successive  increments  of  £,  employing  Equation  (23)  to  extrapolate  to  the  new 
concentration  f(£;+A£)5  7?)  from  previous  f(£j,  tj)  using  forcing  coefficients: 


A,  =  +  (3j) 

Bj  =  (6rtBj-1i/(«iAj.,+fj) 

for  Atj-  =  At}j_j  ,  and 

Aj  =  7/j+i  /  [7j+1  (OjAj.j  +/Jj)  -  7j«j+1  ] 

Bj  =  [7j+1  -7j«j+1]/l7j+i  («jA j.j+Uj)  -7jaj+1] 


(28) 

(29) 

(30) 

(31) 


J+1 


for  At?j  =  l/2Arjj_2 .  The  values  of  fj  are  calculated  from  the  recursion  formula 
i  +  1  i+1 

fj  =  Ajfj+t  +  Bj  (32) 


i+1 


provided  A0,  B0  at  the  centerline  and  fn  at  the  tube  wall  are  both  known. 


b.  Coefficients  at  Centerline,  AQ  and  BQ.  At  the  centerline,  V  ~  d<po/dr}=0 
because  of  symmetry;  diffusion  Equation  (15)  is  reduced  to  the  form 

U3f/3£  =  (I/77SC)  [??32  W  +  0f/3i?] .  (33) 

Equation  (33)  is  indeterminate  at  the  centerline,  as  is  its  corresponding  finite  difference 
form,  because  of  the  indeterminate  terms  (1/rj)  df/drf.  In  the  limit,  as  77-H),  using 
L’Hospital  rule,  this  indeterminate  term  may  be  replaced  by  its  corresponding  determi¬ 
nate  form.  The  centerline  differential  Equation  (33)  thus  has  the  form 

Uaf/a$  =  (2/Sc)(32f/37?2).  (34) 

Noting  the  boundary  condition,  Equation  (17), 

At  7}=0:  df/dy=0  i.c.,f'|1=f,++11.  (35) 
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Equation  (34)  can  be  written  in  the  following  finite  difference  form: 


i+1  i+1 

f()  =  A0^1  +B0  (35A) 

l+'A 

where  A0  =2/  [Sc(Ar?)2  (UQ  /  A£+2/Sc(At?)2  ]  (36) 

i+^2  i  i 

and  B0  =  [U0  /A£-2/Sc(Atj)2  ]  fQ  +  2f1/Sc(Ag2  / 

[l1^2  /  A£+2/Sc(Aij)2  ]  .  (37) 

With  A0  and  B0  calculated  by  Equations  (36)  and  (37)  respectively,  Aj  and  Bj  with 
j  =  1,  2,  —  (n-2)  can  be  evaluated  with  Equations  (28)  and  (29)  or  (30)  and  (31), 
whichever  is  appropriate. 

c.  Wall  Concentration,  fn  .  For  evaluation  of  wall  concentration  fn  ,  we 
begin  with  the  boundary  condition,  Equation  (18): 


Atr?=l:  3f/97j=RemSc(l-f)  . 

Its  corresponding  finite  difference  form  can  be  written  as 


i+1  i+1  i+1 

fn  =fn-l  -AT?ReScm(l-fn.,). 

(38) 

Substituting  Equation  (38)  into  Equation  (23)  with  j  =  n-1  yields 

“n-l  fn-2  +  ^n-1  fn-l  =  5n 

(39) 

where 

(40) 

Cl  =  0,1-1  +(1+AT?ReScm)7n.i 

(41) 

6n-l  =5n-l  +7n-l  AvRcScm. 

(42) 

Elimination  of  fn.2  between  Equations  (39)  and  (32)  with  j  =  n-2  yields  the  following 
solution  for  fn,j : 


i+1 


fn-l  '(5Vl  ‘an-l  Bn-2)/ (an-l  An-2  +^n-l) 


(43) 


Then  f|,+  1 


can  be  calculated  by  Equation  (37);  i.c., 


21 


fn  =  (l+ArjReScm)  fn  l  -  Ar?ReScm,  (44) 

i+1  i+1 

and  fn_j  to  fg  can  be  calculated  by  the  recursion  Equation  (32). 

Thus,  the  problem  is  represented  by  Equation  (43)  for  j=n-l,  Equation 

(44)  for  j=n,  and  Equation  f32)  for  j=n-2,  n-3, - 0,  together  with  initial  condition, 

Equation  (16). 

d.  Sequence  of  Computations  and  Computer  Program.  The  procedure  for 
solution  of  the  diffusion  problem  can  now  be  outlined  by  the  following  sequence  of 
operations: 


(1)  Input  parameters,  Reynolds  number  Re,  Schmidt  number  Sc,  and 
mass  injection  rate  m  are  set. 


(2)  R+  value  corresponding  to  given  Re  is  read  from  Fig.  13. 

(3)  The  entrance  velocity  profile  0(7?)  is  calculated  from  Equation  (9): 

(4)  U,  V,  and  Dg  aic  calculated  by  Equations  (19),  (20),  and  (14), 
respectively. 

(5)  a-,  P-,  7-,  and  6-  are  evaluated  by  Equations  (24),  (25),  (26),  and 
(27)  respectively. 


(6)  A0  and  BQ  are  calculated  by  Equations  (36)  and  (37),  respectively. 

(7)  Aj  and  Bj  with  j  =  1,  2,  —  (»-2),  are  evaluated  by  Equations  (28) 
and  (29)  for  A7jj=Ar?j.j  or  by  Equations  (30)  and  (31)  for  Arjj=l/2Ar?j  j . 

j+1  i+1 

(8)  fn.j  and  fn  are  evaluated  by  Equations  (43)  and  (44). 

i+1 

(9)  fj  ,  with  j  =  (n-2),  (n-3),  — 0,  is  evaluated  by  Equation  (32). 

(10)  The  value  of  i  is  incremented  by  unity,  and  steps  (4)  through  (9) 
are  repeated  until  the  concentration  field  is  computed  for  the  desired  value  of  £. 

A  computer  program  in  Fortran  IV  Language  has  been  written  and  is  ap¬ 
pended  for  the  above  sequence  of  operations.  In  this  program,  Ar^’s  are  specified  in  ac¬ 
cordance  with  the  Table,  and  the  initial  Aj-  was  chosen  to  be  10'6.  Then,  the  successive 
A£  was  twice  of  tb  previous  A£,  or  0.5,  whichever  was  smaller. 
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For  the  convenience  of  the  users  of  the  present  theory,  the  complete 
program  listing,  together  with  users’  instruction,  may  be  found  in  Appendix  B.  Exam¬ 
ples  of  calculated  helium  concentration  field  at  Re  =  1  x  104  to  1  x  106  and  m  =  -0.001 
are  plotted  in  Figs.  15  through  17. 


IV.  RESULTS  AND  DISCUSSION 

9.  Summary  of  Results.  Experimental  and  theoretical  investigations  of  turbu¬ 
lent  diffusion  in  round  tubes  have  been  made. 

Experiments  were  carried  out  with  fully  developed  turbulent  flows  of  air  at 
Re  =  15,000  and  30,000  entering  a  porous  tube  with  circular  cross  section.  Measure¬ 
ments  of  velocity  distribution  were  made  when  air  was  injected  uniformly  through  the 
tube  wall  for  36  radii  at  various  ratios  of  mass  velocity  through  the  tube  wall  to  the 
average  mass  velocity  at  the  entrance  cross  section  of  the  tube  ranging  from  0  and 
0.0004.  These  data  are  plotted  in  Fig.  8. 

Measurements  of  helium  concentration  we;c  also  made  at  Re  =  1.5  x  104  and 
3  x  104  when  helium  was  injected  uniformly  through  the  tube  wall  for  36  radii  at  vari¬ 
ous  ratios  of  helium  mass  velocity  through  the  tube  wall  to  the  mass  velocity  of  air  at 
the  entrance  cross  section  of  the  tube  ranging  from  0.0001  to  0.0004.  These  data  are 
summarized  in  Figs.  9  through  11. 

A  finite  difference  procedure  for  turbulent  diffusion  was  developed  based 
upon  the  mixing-length  concept  for  the  turbulent  diffusivity  and  the  assumption  that  for 
small  wall  injection  the  distortion  of  velocity  profile,  u/0,  due  to  mass  injection  is  negligi¬ 
ble.  A  computer  program  in  Fortran  IV  language  has  been  written  in  accordance  with  the 
developed  calculation  procedure.  Examples  of  calculated  helium  concentration  distribu¬ 
tion  at  Re  =  5  x  103  to  1  x  10s  and  iti  equal  to  -0.001  arc  plotted  in  Figs.  15  through  17. 

However,  the  significance  in  a  finite  difference  solution  of  any  physical  prob¬ 
lem  lies  in  the  accuracy  of  the  finite  difference  scheme  and  the  adequacy  of  the  assump¬ 
tions  made.  The  remainder  of  this  section  is  devoted  to  verifying  the  accuracy  of  the 
present  numerical  scheme  and  physical  assumptions. 

10.  Accuracy  of  Calculation  Scheme.  Accuracy  of  the  calculation  scheme  can  be 
measured  by  comparing  the  solution  with  the  known  exact  solution.  Since  no  exact  so¬ 
lution  is  available  for  turbulent  flow,  tests  for  accuracy  will  first  be  made  by  comparing 
the  present  data  with  the  more  exact  solution  for  laminar  flow23  and  with  correlated 

*°G.  Hailliby,  “Laminar  Heat  Transfer  in  the  Thermal  Entrance  Region  of  Circular  Tubes  and  Two-Dimensional 
Rectangular  Ducts  with  Wall  Suction  and  Injection,”  Int.  J.  Heat  mass  Transfer,  Vol.  14,  pp.  223-243(1971). 
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Fig.  16.  Calculated  helium  concentration  profiles  at  Re  =  5  x  104  and  m  =  -0.001. 
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Helium  concentration 
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Fig.  17.  Calculated  helium  concentration  profiles  at  Re  =  l  x  105  and  in 
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heat  transfer  data  for  zero  wall  injection.24 

The  exact  solution  for  laminar  heat  transfer  with  wall  mass  injection  has  been 
obtained  by  Raithby25  .  In  order  to  compare  his  solutions  with  our  numerical  solution, 
the  present  computer  program  has  been  modified  by  replacing  the  turbulent  velocity 
profile  <Pq  (77)  with  the  well-known  parabolic  velocity  profile  for  laminar  flow  and  re¬ 
placing  the  dimensionless  effective  diffusion,  De,  with  the  laminar  exchange  coefficient, 
I/Pr.  Figure  18  shows  a  comparison  of  the  present  solution  of  temperature  profiles 
with  Raithby ’s  exact  solution  under  the  same  conditions.  Figure  19  shows  a  comparison 
of  Nusselt  numbers  at  the  thermal  entrance  region  calculated  by  the  present  method 
with  those  calculated  by  Raithby. 

1.0 


0.8 


t*  0.6 

I 

\ 

4  °*4 


0.2 


ij(=r/R) 

Fig.  18.  Comparison  of  temperature  profiles  for  laminar  flow  with  constant 
heat  flux  and  uniform  wall  injection:  — , exact  solution  from  Raithby; 

X,  present  theory. 

^W.  M.  Kays,  “Convective  Heat  and  Mass  Transfer,"  McGraw-Hill,  New  York  (1966). 

25 

G.  Raithby ,  “Laminar  Heat  Transfer  in  the  Thermal  Entrance  Region  of  Circular  Tubes  and  Two-Dimensional 
Rectangular  Ducts  with  Wall  Suction  and  Injection,"  Int.  J.  Ileal  Mass  Transfer,  Vol.  14,  pp.  223-243  (1971 ). 
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Fig.  19.  Comparison  of  Nusscll  number  for  laminar  flow  with  constant  heat 
flux:  — ,  exact  solution  from  Raithby;  X,  present  theory. 
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Turbulent  heat-transfer  data  with  zero  wall-mass  transfer  has  been  correlated 
by  Kays.26  In  order  to  compare  the  present  solution  with  Kays’  correlation,  the  present 
program  was  modified  by  replacing  the  Schmidt  number,  Sc,  by  the  Prandtl  number,  Pr. 
Figure  20  shows  a  comparison  of  the  fully  developed  Nusselt  number,  Nu,  calculated  by 
the  present  method  with  that  correlated  by  Kays.  Figure  21  shows  a  comparison  of  the 
Nusselt  number  for  turbulent  heat  transfer  in  the  thermal  entrance  region. 

It  can  be  seen  in  Figs.  18  through  21  that  the  present  results  are  consistently 
in  agreement  with  the  available  laminar  and  turbulent  heat  transfer  data  with  and  with¬ 
out  wall-mass  transfer. 

11.  Justification  of  Main  Assumption  and  Comparison  with  Experiments.  The 
main  assumption  made  in  the  present  analysis  is  that  the  velocity  profile,  u/O,  is  insensi¬ 
ble  to  the  wall-mass  transfer  at  the  range  of  wall-mass  flux,  m,  ranging  from  0.0001  to 
0.0004.  This  assumption  is  justified  by  examining  Fig.  22  where  the  measured  u/u,  de¬ 
rived  from  the  experimental  data  of  Fig.  8  and  the  theoretical  data  of  Equation  (9),  is 
plotted. 


Finally,  the  overall  accuracy  of  the  present  solutions  for  the  turbulent  diffu¬ 
sion  in  round  tube  is  compared  with  the  experimental  data  under  the  same  conditions 
in  Figs.  23, 24,  and  25.  It  can  be  seen  in  these  figures  that  the  agreement  between  the 
experiments  and  theory  is  excellent  with  the  exception  of  the  uncertainy  of  the  Gas 
Analyzer  for  concentration  below  1  percent,  as  pointed  out  by  the  supplier  of  the 
instrument. 


26W.  M.  Kavs,  “Convective  Ilcat  and  Mass  Transfer,”  McGraw-Hill,  New  York  (1966). 
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Fig.  23.  Comparison  of  calculated  and  measured  helium  concentration  at  Re  =  1.5  x  104 
and  m  =  -0.0002:  —  ,  present  theory;  X,  O  and  A,  experiments  at  x/D  =  15,  10,  and  5. 
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V.  CONCLUSIONS 


12.  Conclusions.  The  results  of  this  research  can  je  summarized  as  follows: 

a.  A  method  for  measurement  of  species  concentration  in  a  turbulent 
stream  has  been  developed.  The  method  has  been  applied  to  the  measurement  of  helium 
concentration  in  air. 

b.  A  theory  has  been  developed  for  the  solution  of  diffusion  equations.  A 
simplified  assumption  was  made  for  the  turbulent  velocity  profile.  (The  validity  of  this 
assumj.  ‘ion  for  m  up  to  0.0004  can  be  judged  by  reference  to  Fig.  22.) 

c.  The  theory  is  simple  and  quick  to  use  in  engineering  calculations  because 
the  simultaneous  solution  of  partial  differential  equations  for  conservation  of  mass,  mo¬ 
mentum,  and  species  is  replaced  by  the  solution  of  a  single  partial  differential  equation 
for  conservation  of  species. 

d.  The  accuracy  of  the  present  solution  can  be  judged  by  inspection  of 
Figs.  18  through  21  and  23  through  25,  where  the  present  predictions  are  compared 
with  available  theoretical  predictions  and  the  present  experiments. 

e.  The  present  method  is  restricted  to  the  small  wall-mass  injection.  This 
restriction,  however,  does  not  prevent  its  usefulness  in  a  vast  number  of  engineering 
applications. 

f.  In  the  reverse  osmosis  desalination  of  saline  water,  the  wall-mass  flux  is 
usually  small  and  so  the  present  theory  is  applicable  with  slight  modification  of  bound¬ 
ary  condition  (Equations  (2)  through  (5)).  Thus,  the  present  theory  should  be  used  for 
optimization  of  the  reverse  osmosis  process. 

g.  In  the  use  of  a  catalytic  converter  for  control  of  pollutant  emission,  the 
net  wall-mass  flux  is  zero  and  so  present  theory  is  again  applicable.  However,  the  reac¬ 
tion  rate  at  the  converter  surface  will  determine  the  boundary  condition  for  this  problem. 
The  present  theory  can  be  used  to  predict  diffusion  of  pollutants  in  the  gas  stream,  and 
subsequently  performance  of  the  converter  can  be  calculated. 

h.  In  application  to  the  water  supply  system  where  dissolution  of  minerals 
and  its  salt  must  be  predicted,  the  present  theory  should  be  used  to  calculate  the  water 
quality. 

i.  In  some  engineering  devices,  for  example  gas  turbine  combustors  and  d/f- 
fusers,  the  wall-mass  injection  is  large.  Under  such  circumstances,  the  present  theory 
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should  be  extended  to  the  regjon  of  large  wall-mass  injection.  The  conservation  equa 
tions  for  continuity  of  mass,  momentum,  and  species  will  have  to  be  solved  simultane 
ously.  In  view  of  the  fact  that  the  momentum  and  diffusion  equations  are  both  para¬ 
bolic,  the  present  method  can  readily  be  extended  to  the  problems  with  large  wall-ma 
injection. 
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APPENDIX  A 


ENTRANCE  VELOCITY  PROFILE,  <p0  (r?) 


To  derive  the  expression  for  velocity  profile  at  the  entrance  to  a  porous  tube, 
4>0(v),  analyses  begin  with  the  differential  equation 


-pc  (du/dr)  /  Tw  =  r/r,v  =  7? .  (Al) 

In  this  equation,  the  effective  viscosity  for  turbulent  flow,  pc,  is  defined  as  the  sum  of 
the  laminar  and  turbulent  viscosities,  p  and  pt: 


MC=M  +  Mt. 


(A2) 


Among  the  various  candidate  representations  for  the  turbulent  viscosity ,  a  model  involv¬ 
ing  Prandtl’s  mixing  length  and  Van  Driest’s  damping  factor  has  gained  broad  acceptance 
in  recent  years.27  '  32 

This  approach  is  employed  here  and  gives  the  following  model  for  the  turbulent 
viscosity: 


pt  =  pi2  /3u/3r/ 

(A3) 

1/R  =  DF  (0.14  -  0.08j?2  -  0.06tj4) 

(A4) 

DF  =  1-cxp  [  -R \/tw/p  (l-i?)/Af]  . 

(A5) 

27It  B.  Kinncv  and  E.  M.  Sparrow,  “Turbulent  Flow,  Heat  Transfer,  and  Mass  Transfer  .n  a  Tube  with  Surface 
Suction,”  ASME  Paper  No.  69-1  IT-4  (1969). 

28E.  M.  Sparrow,  V.  K.  Jonsson,  G.  S.  Beavers,  and  It.  G.  Over,  “I ncomprcssiblc  T ur)  ’t  Flow  in  a  Permeable 

Wall  Duct,”  A$ME  Paper  No.  71-WA/FE-J  (1971). 

29T.  Ccbcci,  A.  M.  0.  Smith,  and  G.  Mosinskis,  “Solution  of  the  Incompressible  Turbulent  Boundary-Layer 
Equations,"/.  Heat  Transfer,  Vol.  92,  pp.  133-143  (1970). 

30S.  V.  Patankar  and  D.  B.  Spalding,  “Heat  and  Mass  Transfer  in  Boundary  Layers,”  Morgan  and  Grampin, 
London  (1967). 

31S  \v  Chi  and  C.  C.  Cliang,  “Effective  Viscosity  in  a  Turbulent  Boundary  Layer,”  AIAA  Journal,  Vol.  7,  pp. 
2032-2035(1969). 

32W.  J.  Glowackr  and  S.  W.  Chi,  “Effect  of  Pressure  Gradient  on  Mixing  Length  for  Equilibrium  Turbulent 
Boundary  Layers,"  AIAA  paper  No.  72-213(1972). 
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Equation  (A3)  expresses  Prandtl’s  mixing  length  concept,  while  Equation  (A4)  is  a 
representation  for  the  mixing  length,  1,  deduced  by  N’ikuradse33  from  measured  velocity 
profiles  for  the  circular  tube.  Van  Driest’s  damping  factor34  for  an  impermeable  wall  is 
expressed  by  Equation  (A5).  The  quantity  A  is  the  damping  constant  which  is  equal  to 
26. 


The  effective  viscosity  expressions  of  Equations  (A2)  through  (A5)  may  now  be 
substituted  into  Equation  (Al)  yielding  the  following  governing  equation  for  the  non- 
dimensional  axial  velocity,  u+  (=  puy/rlp/ p): 


t?  =  -  (1/R+)  (du+/d7?)  +  {  (0.14-0.08t?2-0.06i?4) 


[1-exp  (-R+(l-i?)/A)]  [  (du+/dr})2. 


(A6) 


Integration  of  Equation  (A6)  with  the  boundary  condition  u+(R+)  =  0  followed  by 
considerable  rearrangement  gives  the  following  expression  for  the  axial  velocity 
distribution: 

u+  =  I"1  2R+tj  {  1  +  [  1+4tj(R+)2  {  (0.14-0.08t?2-0.06t?4) 

(1-exp  (-R+  (l-7?)/A)}2  }Vl  ]•  ‘dr? .  (A7) 

In  addition,  the  following  relations  can  be  verified  in  a  straightforward  manner: 

Rq  =  2  j^Tjr/dr?  (A8) 

Re  =  2R+Uq  (A9) 

0O  u+/u^  (A10) 

Using  Equations  (A7)  through  (A9),  <Pq(v)  at  different  Re  values  and  Re  versus  R+ 
relation  can  be  evaluated. 


33tf.  Schlichtiiw,  “Boundary  Layer  Theory,”  McGraw-Hill,  New  York  (1968). 

^E.  R.  Van  Driest,  "On  Turbulent  Flow  Near  a  Wall,”/.  Aero.  ScL,  Vol.  23,  pp.  1006-1011, 1036(1956). 


y 


39 


APPENDIX  B 


LISTING  AND  USERS’  INSTRUCTION  OF  COMPUTER  PROGRAM 
FOR  TURBULENT  DIFFUSION  IN  ROUND  TUBES 


A  computer  program  for  calculation  of  turbulent  diffusion  in  round  tubes  has 
been  written  in  Fortran  IV  language.  The  program  has  been  checked  out  on  CDC  6000 
series  and  IBM  360  series  computers.  A  complete  listing  of  the  computer  program  is 
presented  below,  and  main  FORTRAN  input  and  output  symbols  are  defined.  All  the 
input  integers,  beginning  with  L  through  N,  are  5-digit  figures  (15),  and  all  other  input 
constant  are  15-digit  figures  (E15.6). 

Input  Symbols: 


» 


TSTA 

a  constant  to  suppress  underflow  1  x  10"2S 

TSTF 

a  constant  to  suppress  underflow  1  x  10‘2S 

NRPS 

number  of  nondimensional  radius  R+ 

RPS 

nondimensional  radius,  R+ 

SC 

Schmidt  number 

SCT 

turbulent  Schmidt  number 

NZETA 

number  of  increment  in  axial  direction 

DDZT 

initial  increment  in  axial  direction,  5jj,  1  x  10'6 

NMBAR 

number  of  nondimensional  wall-mass  injection  rate 

AMBAR 


nondimensional  wall-mass  injection  rate,  m 


Output  Symbols: 


RE  Reynolds  number  at  entrance  to  diffusion  section 

ETA  nondimensional  radius,  i\ 

DU/DETA  nondimensional  gradient  of  velocity  d0Q/dr? 

U/UB  nondimensional  axial  velocity  defined  as  local  velocity  divided  by  bulk 

average  velocity,  <j>0 

U/UC  nondimensional  axiai  velocity  defined  as  local  velocity  divided  by  center- 

line  velocity 

SC  Schmidt  number 

SCT  turbulent  Schmidt  number 

MB  nondimensional  wall-mass  injection  rate,  m 

ZETA  nondimensional  axial  position,  £  (•-  x/R) 

F  mass  concentration  of  species  under  consideration,  f 

FBAR  species  concentration  defined  as  (fw--f)/(fw-f) 

FCB  species  concentration  defined  as  (f-f^)/ (f w-f0) 

Nu  Nusselt  number 
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PROGRAM  MAIM  ( INPUT , OUTPUT , TAPE5=INPUT , TAPE6=0UTPUT) 

C  TURBULENT  OIFFUSION  IN  ROUND  TUBES 

DIMENSION  ETA  (241)  ,  U0(241),  V(241),  F(241i,  A(241),  8(241'), 
1U1  (241 )  ,  U2  (241)  ,  DUO  (241)  ,  EDEK241),  EDE2(241) 

WRITE  (6,21) 

READ  (5,12)  TSTA,  TSTF 
READ  (5,1)  NRPS 
00  101  NR=1, MRPS 
REAO  (5,11)  RPS 
ETA ( 1) =0  . 

UO(1>=0. 

OUO(l)  =0. 


1 

FORMAT 

(15) 

11 

FORMAT 

(E15.6) 

12 

FORMAT 

(2E15.6) 

13 

FORMAT 

(3E15.6) 

14 

FORHAT 

(4E15.6) 

21 

FORMAT 

(//35H  TURBULENT 

OIFFUSION  IN  ROUND 

TUBES) 

22 

FORMAT 

( /4H  RE= , E15. 6) 

23 

FORHAT 

( /53rl  ETA 

OU/OETA 

U/UB 

U/UC) 

24 

FORMAT 

(/53H  RE 

SC 

SCT 

MB) 

25 

FORMAT 

(/6H  ZETA= ,E15 .6) 

26 

FORNAT 

{/53H  ETA 

F 

FBAR 

FCB) 

27 

FORMAT 

(/5H  NU= , E 15 . 6) 

1=1 
J  =  0 
N=1 

DO  102  IM=1, 8 
IF  (N-l)  211,201,211 

211  IF  (N-2)  Z12 ,202,212 

212  IF  (N-3)  213,203,213 

213  IF  (N-4)  214,204,214 

214  IF  (N-5)  215,205,215 

215  IF  (N-6)  216,206,216 

216  IF  (N-7)  217,207,217 

217  IF  ( N-8)  218,208,218 

218  CONTINUE 

201  DET  A  =  0 . 00  64 
M=4  0 

GO  TO  220 

202  DET A  =  0 .0  032 
M=40 

GO  TO  220 

203  DETA=0 .0016 
M=40 

GO  TO  220 

204  OET A=0 . 0008 
M=40 

GO  TO  220 

205  OETA=Q . 0004 
M=40 

GO  TO  220 

206  OETA=0 .0002 
H=10 

GO  TO  220 

207  OETA=0 . 0001 
M=10 

GO  TO  220 

208  OET A=0 , 00  00 •> 

M=20 

GO  TO  220 
220  00  103  K=i,M 
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J=J+1 

1=1+1 

IF  (K-l)  222,223,22; 

223  IF  (N-l)  222,224,22 

224  E1=0 • 0064 
E2=0 . 0128 
ETA (I) =0 . 0128 
GO  TO  225 

222  E1=ETA ( J) +OETA 
E2=E1+0ETA 
ET  A  ( I)  =  E2 

225  0U1=2. *RPS* El/ ( l.+SQRTU. +4. *E1*RPS*RP3* ( 0.14-0. 080 +E1++2- 0.06 
i*El**4)  +*2*(1. -EXP  ((-l.+RPSM  1-El) ) /26. )  >**2)> 

0U2=2.*RPS*E2/ <1.+SQRT(1.+4.*E2*RPS*RPS* ( 0 . 14-0 . 060+E2++2-0. 06 
l*E2+*4)**2*(l.-EXP((-i.+RPS*(l-E2) )/26.) )+*2)) 

DUO (I) =OU2 

UO(I)=DETA*(DUO(J> +4. +0U1+0U2) /3. 

UO(I)=UO(J)+UO(I) 

103  CONTINUE 
N=N+1 

102  CONTINUE 

00  104  1=1,24 
K=I-1 

UO ( I) =UO (241)— UO(I> 

J=I  +  1 

IF  (1-1)  105,106,105 
106-  U8=0. 

GO  TO  104 
105  CONTINUE 

107  UB=UB+ (ETA(I)*UO(I)+ETA(K)*UO(K)  !  * (ETA { I ) -ETA t K) ) 

104  CONTINUE 
RE=UB*RPS 
R£H=2.*RE 
J=9 

WRITE  (6,22)  REW 
HRITE  (6,23) 

00  108  1=1,241 
J=J  +  1 

UO ( I) =UO  < I) /UB 

•  OUO(I)=OUO(I)/UB 

•  UUC=UO(I)/UO(l) 

IF  (J-10)  138,110,108 
110  HRITE  (6,14)  ETA(I) , OUO (I ) , UO ( I ) ,UUC 
J=0 

108  CONTINUE 

REAO  (5,13)  SC,SCT 
REAO  (5,1)  NZETA 
REAO  (5,11)  OOZT 
REAO  (5,1)  NMBAR 
ZET  A  =  0 

00  500  IM?=1,NMPAR 
REAO  (5,11)  AMBAR 
WRITE  (6,24) 

REW=2. *DE 

WRITE  (6,14)REW,SC,SCT,AMPAR 
ZETA=0. 

OZET  A=OPZT 
DO  450  INZ=1, NZETA 
IF  (INZ-23)  871,871,872 
371  DZET  A=2 .  +  OZETA 
GO  TO  873 

872  OZET A  =  0. 5 

873  CONTINUF 

IF  (INZ-1)  451,452,451 
452  CONTINUE 

00  440  K=l, 241 
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F(K) =0 

IF  (K-i)  456,457,456 
457  V(K) =0 • 

V8=0. 

VVK=0. 

GO  TO  611 
456  KA=K-1 

VC=RE*AM8AR*U0(K)*ETA(K) 

VVK  =VVK  +1VB+VC) *(ETA(K)-ETA(KA)) 

V (K) =VVK/£TA (K) 

V  B=VC 

511  CONTINUE 

Ui <K) =UO (K) *RE 
AKI=AMBAR*R£*(l.-ETA(K))/2 

0F=SQRT(AKI**4+4.* (RPS*( 1 .-ETA <K) )/26.)**4> 
OF=$QRT(OF+AKI**2)/(SQRT<2.) > 

OF=l-EXP(-AKI-OF) 

IF  (K-i)  613,612,613 

613  CONTINUE 

EOE1 CK)=OF* (0. 14-0.08*ETA(K) **2-0. 06*ET A < K> **4> 

E0E1(K)=RE*0U0(K)*E0E1(K)**2/SCT 

EOE1 (K) =1 . /SC+EOE1 (K) 

EOE1 (K) =ET  A ( K) *£0E1(K) 

GO  TO  440 
612  EOE1  (K) =0 • 

440  CONTINUE 
GO  TO  454 

451  DO  453  K=l, 241 
U1(K) =U2 (K) 

EOEi (K) =EDF2 (K) 

453  CONTINUE 

454  ZETA1  =  Z£-  A 

ZET A=ZET  i+OZE f  A 
DO  455  K-i, 241 

U2(K)=RE*(1.-2.*AHBAR*ZETA)*U0(K) 

IF  ( K— 1 )  614,615,614 

614  CONTINUE 

E0E2(K)=E0El(.O/ETA(K)-l./SC 

EDE2(K)=EOE2(K)* < 1 .-2 . *AMBAR*ZET A ) /( 1 .-2 . * AHBAR* ZET Al> 
EOE2<K)  =  (EDE2(K) *1 ./SC) *  ETA ( K) 

GO  TO  455 

615  EOE2 (K)  =  0 . 

455  CONTINUE 

DO  301  IN=i , 6 

IF ( IN- 1)  411,401,411 

411  IF ( IN- 2)  412,402,412 

412  IF( IN-3)  413,403,413 

413  IF ( IN-4)  414,404,414 

414  IF ( IN-5)  415,405,415 

415  IF (IN-6)  416,406,416 

416  IF (IN-7)  417,407,417 

417  IF (IN-8)  418,408,418 

418  CONTINUE 

401  OET  A  =  fl ,0128 
M=4  0 

GO  TO  420 

402  OETA=0 .0064 
H=40 

GO  TO  420 

403  OET A  =  0 .00  32 
M=40 

GO  TO  420 

404  DET A=0 . 0016 
M=4  0 

GO  TO  420 

405  OETA=0 . 0008 
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m=4u 

GO  TO  429 
406  OETA=O.ODC4 
M=10 

GO  TO  420 
40/  OETA=0. 0902 
H=10 

GO  TO  420 
400  OET A  =  0 .0001 
H=20 

420  IF  (IN-i)  431.430,431 

430  PET A  =  0 .5* (U1 ( 1 ) +  U2 ( 1 ) ) /DZETA+2./ ( SC* DETA+DETA) 
GAHA=-2./(SC*0ETA»DETA) 

DELT=0.5+ (U1 ( 1 ) +U2 (1 ) ) 'OZETA-2.7 (SC*OETA»OETA> 
DELT=OELT*F(l) + 2 . *F (2 ) / , SC*CETA*9ETA ) 

A(l) =-GAHA/PET  A 
B(l) =OELT/?ETA 
1  =  3 

MMzQ 

431  00  441  K=l,3 
1  =  1  +  1 

J  =  I  +  1 
L=J  +  1 

IF  K-4)  433,436,433 

436  IF  (IN-8) 437,434,437 

437  L=L  +  1 

433  E01=9.25»(C:E1(I)+E0E1(J)+E!)E2<I)+EDE2(J) ) 

E02  =  9*  26+  (E'JEl  <J>+E0E1(L)+'0E2(J>+EDE2(L> ) 
UU=0.5»(U1(J)+U2(J)) 

VV=0.5*(V(J)+\MJ)) 

ALFA=-n.5*E91/ (ETA(J)+OETA*O£TA)-0.25*VV/OETA 
8ETA=UU/OZETA+9.5* (E01+E02) / (ETA (J) *D£TA*OETA) 
GAMA=0.25*VV/OETA-0.5*FO2/(ETA(J)*OFTA*OETA) 
OELT=F(I)*(0.5*E91/(£TA<J)+O£TA*OETA>+9.25*VV/DETA) 
l  +  F(J)  +  {U'J/0ZETA-a.5ME01  +  E02)/(ETA(J)*0ETA+0ETA)  ) 
2+F(L)*(P.5+EO2/(ETA(J)»OETA*OETA>-0.25*VV/DETA) 

IF  (L-J-l)  439,625,438 
625  IF  (HM)  470,439,470 
439  A ( J) =-GAMA  / ( A LF A  * A ( I ) +0ET A ) 

B(J>=(0ELT-ALFA*9(I) ) /(ALFA* A (I) +PETA) 

GO  TO  441 

438  I=J 
J=J  +  1 
HM=31 

OETA=0.5*DETA 
ALF  Al  =  ALC A 
8ETA1=6ET  A 
GAMA1=GAHA 
DEL  T l  =  OEL  T 
GO  TO  433 
470  MM=I-1 

A(I)=GAMA1*9ETA/(GAMA+ (ALFA1+ A (MM) +FETA1) -GAHA1* ALFA) 
8(1)  =  (GAMAM0ELT1-ALFA1*3(MM)  )  -G  AKA  1  *DEL  T)  /( GAMA* 

1  (ALFA1  +  A (MM) +8ETA1 ) -G A“A1 +ALFA) 

I=MM 
MM=  0 

GO  TO  441 

434  AA4=0ETA*PE*SC*AH9A9 
BETA=9ETA+GAMA/(1.-AA«) 

OEUT=OElT  +  AAM*GAMA/(l.-AA“i) 

11=1-1 

F(I)  =  (QFLT- ALFA*  Kill  ) / ( ALFA *6 ( 1 1) +rETA ) 

F (J)  =  (F (I)  —  A4M)/ (1,-AAW) 

441  CONT  l'lj* 

301  CONT IMUr 

00  -,74  <  =  !,■>  n 
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I  =  2  4  0  -  K 
J=l  +  1 

IF  (A-S(A(I))-TSTA>  731,792,702 

702  IF  (A>)S(P(J)  ) -TSTC)  7?i,  7'T,  703 

703  F(I>=ACI> *F(J) +0(1) 

GO  TO  (,7o 

701  F<I)=0(I> 

470  CONTINUE 

WRITE  (6,25)  ZETA 
FR=F(24i) 

FB=-2**ZETA  +  AM)AR/'{1,-2.*Z£7A*AM!)AR) 

RNU  = (FR-F (240) )/ (0ETA+ (FR-FO) ) 

RNU=2.*RNU 

WRITE  (6,27)  RNU 

WRITE  (6,26) 

00  475  1=1,241,10 
FBAR=(FR-F(I) ) /(FR-PR) 
FC8=(F(I>-F(241)  )/(F(l)-F(241)) 

WRITE  (6,14)  ETA(I),F(I),F8AR,FCB 
475  CONTINUE 
450  CONTINUE 
02ETA=002T 
500  CONTINUE 
101  CONTINUE 
CALL  EXIT 
STOP 
END 
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36/26/7 
14.  in .  (.6, 

14.10.46. 

14.10.47. 

14.10 .47. 

14.10.48. 

14.10.48. 

14.10.49. 
14.10.49. 
14.10.49. 


JtapST  3*3P1  M£R0C 

MTAP3.CM60000.T400.  PATTI 
TASK(TNMT72039,PW****»,TRTS  <jrTO> 

COPYCF (INPUT, X)  ''RIS.SCIR) 

REHINO(X) 

COPYSOF(X)OUTPIJT) 

CP  003.033  sec. 

PP  001.390  sec. 

ooa.081  sec. 
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